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Abstract 

This  paper  examines  a  buffer  scheme  to  mitigate  the  negative  impacts  of  power-conditioned  loads  on  network  voltage  and  transient  stabilities. 
The  scheme  is  based  on  the  use  of  battery  energy-storage  systems  in  the  buffers.  The  storage  systems  ensure  that  protected  loads  downstream  of 
the  buffers  can  ride  through  upstream  voltage  sags  and  swells.  Also,  by  controlling  the  buffers  to  operate  in  either  constant  impedance  or  constant 
power  modes,  power  is  absorbed  or  injected  by  the  storage  systems.  The  scheme  thereby  regulates  the  rotor-angle  deviations  of  generators  and 
enhances  network  transient  stability.  A  computational  method  is  described  in  which  the  capacity  of  the  storage  systems  is  determined  to  achieve 
simultaneously  the  above  dual  objectives  of  load  ride-through  and  stability  enhancement.  The  efficacy  of  the  resulting  scheme  is  demonstrated 
through  numerical  examples. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  stability  of  power  systems  has  always  been  an  important 
consideration,  especially  when  it  involves  the  transmission  of 
power  over  long  distances.  The  problem  has  been  compounded 
somewhat  in  recent  years  as  a  significant  proportion  of  present- 
day  industrial  loads  are  equipped  with  power-conditioning 
capabilities.  These  loads  usually  draw  constant  power  regard¬ 
less  of  their  terminal  voltage  level.  In  the  event  of  a  fault  in  the 
power  system,  the  depressed  voltage  will  result  in  an  increase  in 
the  load  current.  This  would  lead  to  an  increase  in  the  voltage 
drop  across  the  source  impedance,  which  in  turn  causes  a  further 
decrease  in  load  terminal  voltage.  Load  voltage  stability  would 
be  impacted  negatively  if  no  countermeasures  were  taken.  Fur¬ 
thermore,  constant-power  loads  also  affect  rotor-angle  stability 
of  power  systems,  in  a  manner  that  can  be  very  much  different  if 
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the  loads  were  to  behave  like  a  constant-impedance  type.  Hence, 
it  is  not  surprising  that  load  modelling  and  the  impacts  of  loads 
on  power  systems  are  attracting  the  interest  of  many  researchers 

[1-3]. 

In  an  attempt  to  address  the  above-mentioned  stability  prob¬ 
lem  due  to  the  presence  of  constant-power  loads,  the  authors  have 
proposed  [4]  the  use  of  a  power  buffer  system  through  which  the 
loads  are  supplied.  The  buffer  consists  of  two  back-to-back  con¬ 
verters  interconnected  by  a  DC-link.  During  a  fault,  the  input 
impedance  (viewed  from  the  source  side)  of  the  buffer-load 
combination  is  controlled  by  the  front  converter  to  maintain 
an  impedance  value  identical  to  that  under  the  pre-fault  con¬ 
dition.  In  this  way,  the  threat  of  an  increase  in  the  load  current 
during  the  fault  event  is  alleviated.  An  energy- storage  system 
was  connected  across  the  DC-link  [4]  and  its  role  was  to  sup¬ 
ply  the  shortfall  in  real  power  between  that  supplied  from  the 
upstream  system  and  the  (constant)  load  demand.  The  buffer 
design  was  extended  [5]  to  cater  for  unbalanced  faults,  with  the 
buffer  controller  designed  to  mitigate  negative  phase  sequence 
voltage  under  the  fault  condition.  The  scheme  allows  not  only 
the  load  to  ride  through  the  fault  event  but  also  alleviates  the 
threat  of  voltage  instability. 
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The  success  of  such  a  buffer  system  must  include  a  suitable 
energy- storage  system  to  compensate  for  any  of  the  mismatch  in 
power  between  the  constant-load  demand  and  that  supplied  by 
the  upstream  system,  during  and  following  the  fault  disturbance. 
This  aspect  of  buffer  design  was  not  analyzed  in  reference  [4]. 
Indeed,  the  energy- storage  capacity  aspect  of  the  buffer  scheme 
is  one  main  design  factures  to  be  addressed  in  the  study  pre¬ 
sented  here.  The  design  will  be  expanded  to  a  proposed  network 
stabilization  scheme.  It  will  be  shown  that  the  scheme  is  capable 
of  extending  the  function  of  the  buffers  described  in  references 
[4,5];  it  can  help  in  achieving  the  dual  purpose  of  network  tran¬ 
sient  stability  enhancement  while  maintaining  supply  quality  to 
the  load.  The  examined  power  system  is  described  in  Section  2. 
The  operational  principle  of  the  buffer  scheme  and  a  method  to 
determine  the  capacity  of  the  storage  system  necessary  to  ensure 
the  success  of  the  scheme  are  also  given.  Numerical  examples 
are  used  in  Section  3  to  illustrate  the  application  of  the  proposed 
scheme. 

2.  Model  development 

Although  practical  power  systems  are  highly  complex,  the 
generic  two-generator  power  system  network  shown  in  Fig.  1 
is  convenient  to  use  for  illustrating  the  operating  principle  of 
the  buffer  scheme.  Indeed,  there  are  radial  networks  that  can  be 
approximated  by  the  network  model,  as  shown  in  reference  [6]. 
One  can  consider  the  practical  situation  in  which  a  group  of  gen¬ 
erators,  represented  by  the  equivalent  unit  G2,  is  interconnected 
with  another  group  of  generators,  represented  by  an  equivalent 
unit  G\ ,  through  the  transmission  link  between  Buses  1  and  2. 
The  power  outputs  of  G\  and  G2  are  denoted  as  Pgx  and  Pq2, 


respectively.  Losses  in  the  transmission  system  are  assumed  to 
be  negligible.  As  the  focus  of  this  work  is  on  the  buffer  scheme, 
details  of  the  generators  and  their  control  systems  will  not  be 
described  here;  although,  in  Section  3  the  illustrative  examples 
do  include  detailed  model  representations. 

In  the  present  problem  formulation  where  power  electronics- 
based  devices  have  been  increasingly  used  in  the  power  systems, 
one  recognizes  some  load  demands  would  be  considered  essen¬ 
tial.  The  nature  of  these  demands  would  be  such  that  supply 
to  them  must  not  be  disrupted  under  specified  momentary  volt¬ 
age  disturbance  conditions.  Indeed,  consideration  for  supplying 
such  loads  has  given  rise  to  the  idea  of  Custom  Power  and  the 
concept  of  a  Power  Quality  Control  Centre  [7,8].  In  whatever 
form,  provisions  would  have  to  be  made  in  the  networks  so  that 
these  loads,  denoted  as  S \  and  S2  in  Fig.  1,  can  ride-through 
such  voltage  sag-swell  disturbances.  Hence,  Si  and  S2  would 
draw  constant  real  powers  Py,  and  Ps2,  respectively,  regard¬ 
less  of  the  status  of  the  power  system.  To  ensure  Si  and  S2  real 
power  demands  are  met,  the  provision  considered  in  this  work  is 
through  the  arrangement  of  the  buffers,  as  shown  in  the  figure. 
More  will  be  said  about  the  buffer  scheme  later.  The  remaining 
loads  D\  and  £>2  at  the  respective  Buses  1  and  2  have  no  such 
stringent  power  requirement.  The  steady- state  powers  drawn 
by  these  loads  are  denoted  as  Ppl  and  Pp2 •  T>1  and  D2  could 
be  assumed  to  be  of  constant-impedance  type  in  this  network 
model. 

Having  met  the  demands  Psl  and  Ppl ,  it  is  assumed  Bus  1 
has  surplus  generation  from  Gi ,  while  Bus  2  has  insufficient  in- 
feed  from  G2  to  meet  the  demands  Ps2  and  Pd2  •  The  deficiency 
is  imported  from  Gi ,  via  the  link,  by  the  amount  P\2.  Therefore, 
the  steady-state  power  flow  pattern  is  as  shown  in  Fig.  1. 


vG2ze2 


Fig.  1.  Schematic  diagram  of  two-generator  network  system  with  power  buffers. 
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2.1.  Buffer  scheme  to  achieve  improved  load  ride-through 

Next,  consider  the  function  of  the  buffers.  In  the  proposed 
scheme,  the  power  buffers  can  be  based  on  electronically  con¬ 
trolled  power  converters.  Under  healthy  system  conditions,  the 
buffers  operate  under  the  constant-power  mode  to  meet  the 
demands  Psl  and  Ps2,  through  control  of  the  respective  front- 
end  converter  while  maintaining  the  DC-link  voltages,  r>DCi  and 
r>DC2>  at  their  nominal  levels.  This  will  ensure  proper  opera¬ 
tion  of  the  downstream  converter-loads  Si  and  £2.  It  is  through 
the  control  actions  of  the  converters  that  Si  and  S2  will  enjoy 
excellent  supply  quality — an  attribute  similar  to  that  described 
in  reference  [8]  for  so-called  super-premium  loads. 

When  a  fault  occurs  in  the  transmission  link,  the  voltage 
in  the  network  will  be  depressed.  To  avoid  the  negative  influ¬ 
ence  of  constant  power  loads  on  system  voltage  stability,  the 
power  buffers  may  be  switched  into  constant  impedance  mode 
by  adjusting  the  PWM  control  signal  [5] .  Under  such  a  condition, 
however,  the  power  supplied  from  the  utility  to  the  downstream 
constant  power  loads  will  be  reduced.  Therefore,  the  impacts  of 
the  constant  demands  Psl  and  Ps2  on  the  buffers  should  also  be 
considered. 

As  stated  above,  over  these  instances  when  the  buffer(s)  have 
to  operate  under  constant-impedance  mode  while  demands  Psx 
(or  Ps2)  have  to  be  satisfied,  the  shortfall  between  P$x  (or  Pg2) 
and  the  power  supplied  by  the  upstream  utility  has  to  be  met.  This 
is  achieved  by  drawing  power  from  the  energy-storage  systems 
within  the  buffer  systems.  Several  possible  storage  media  such 
as  supercapacitors  and  flywheels  could  be  used.  In  this  investi¬ 
gation  a,  lead-acid  battery  energy- storage  system  (BESS)  has 
been  considered.  This  is  because  the  lead-acid  battery  is  one 
of  the  most  cost-effective  and  well-established  technologies  for 
large-scale  energy- storage  [9].  The  BESSs  are  shown  connected 
across  the  DC-link  of  the  buffers  in  Fig.  1. 

2.2.  Buffer  scheme  to  achieve  stability  enhancement 

With  the  above-described  power  buffers  installed,  their  func¬ 
tion  can  be  extended  to  act  as  load  dampers  or  stabilizers  [5]. 
The  operation  scheme  of  the  power  buffer  in  enhancing  system 
transient  stability  is  addressed  as  follows. 

From  the  instant  of  the  fault  inception  and  the  initial  transient 
period  of  (say)  a  few  seconds,  the  governor  control  systems 
of  G\  and  G2  will  generally  be  too  slow  to  respond  appre¬ 
ciably.  Thus,  the  input  mechanical  powers  to  the  generators 
can  be  considered  constant.  On  the  other  hand,  the  excitation 
systems  would  usually  be  able  to  adjust  the  respective  gen¬ 
erator  internal  emf,  since  the  control  systems  are  designed  to 
regulate  the  terminal  voltage  of  the  generators.  Even  with  the 
actions  of  the  excitation  systems,  however,  a  severe  fault  would 
cause  the  power  P\2  transmitted  through  the  link  to  be  greatly 
reduced.  If  both  Si  and  D\  were  to  behave  like  a  constant- 
impedance  type,  Bus  1  load  demand  would  also  reduce  as  the 
load  terminal  voltage  is  depressed  during  a  fault.  Therefore, 
the  surplus  energy  resulting  from  the  constant  input  mechan¬ 
ical  power  and  the  reduced  output  electrical  power  Pg1  will 
accelerate  the  rotor  of  Gi .  To  reduce  the  power  imbalance  at 


Bus  1  and  since  Si  is  of  a  constant-power  type,  it  is  proposed 
that  when  the  Gi  rotor  accelerates  during  a  fault,  Power  Buffer 
1  will  be  designed  to  operate  under  a  constant-power  mode. 
In  this  way,  the  load  demand  at  Bus  1  will  be  higher  than 
that  without  the  buffer.  This  will  reduce  the  acceleration  of 
G\. 

To  change  the  buffer  operation  mode  based  on  Gi  rotor  speed 
may  pose  some  implementation  difficulty  since  Gi  is  an  equiv¬ 
alent  of  a  cluster  of  generators.  The  precise  meaning  of  rotor 
speed  will  need  further  careful  consideration.  As  an  alternative 
to  rotor  speed,  the  transfer  power  P\2  can  be  used  instead  to 
identify  the  case  of  acceleration  or  deceleration  of  Gi .  When  a 
fault  occurs  and  P\2  reduces,  Gi  accelerates  and  the  operation 
mode  of  Buffer  1  is  switched  into  constant-power  mode  to  reduce 
the  acceleration.  Conversely,  when  P\2  increases  above  its  pre¬ 
fault  level,  Buffer  1  is  then  switched  into  constant-impedance 
mode. 

At  Bus  2,  the  situation  is  the  reverse  of  that  of  Bus  1 :  when 
P\2  reduces  and  without  Buffer  2,  G2  would  decelerate  with 
respect  to  G\ .  This  is  the  classical  2-machine  transient  stabil¬ 
ity  situation.  Depending  on  the  pre-fault  power  transfer  level 
and  fault  clearing  time  (among  other  factors),  there  is  a  dan¬ 
ger  Gi  and  G2  may  go  into  a  first-swing  instability  condition. 
Hence,  to  mitigate  the  possibility  of  instability,  Buffer  2  will 
be  designed  to  operate  under  a  constant-impedance  mode  dur¬ 
ing  a  fault  so  that  Bus  2  total  load  demand  will  be  reduced. 
In  this  way,  the  deceleration  in  G2  would  be  less  compared 
with  the  case  without  Buffer  2,  because  in  the  latter  case,  the 
constant-power  demand  Ps2  has  to  be  met  by  an  increase  in 
Pq2  since  P\2  is  reduced.  Thus,  by  forcing  Buffer  1  into  a 
constant-power  mode  and  Buffer  2  into  constant-impedance 
mode  during  a  fault,  the  rotor  speed  difference,  and  there¬ 
fore  the  rotor-angle  difference,  between  G\  and  G2  can  be 
reduced  compared  with  to  the  case  when  the  buffers  are  not 
used. 

The  buffers  can  also  be  designed  to  operate  on  the  same  prin¬ 
ciple  following  fault  clearance  when  the  network  readjusts  its 
bus  voltages  and  power  flows.  When  P\2  reduces  and  there  is 
accelerating  power  acting  on  G\  to  increase  the  generator  speed, 
Power  Buffer  1  is  to  operate  under  a  constant-power  mode.  Con¬ 
versely,  when  P\2  increases  and  the  generator  decelerates,  the 
buffer  is  switched  to  a  constant-impedance  mode.  The  role  of 
Buffer  2  is  to  regulate  the  speed  of  G2  in  a  similar  manner. 
Repeated  switching  between  the  two  buffer  modes  continue  until 
the  network  settles  into  a  new  steady- state  operating  condition. 
Hence,  it  would  appear  that  P\2  variations  can  be  used  to  control 
effectively  the  operating  modes  of  the  buffers. 

2.3.  DC-link  voltage  constraint 

In  the  scheme,  the  buffers  function  to  guarantee  a  constant 
power  supply  to  £j  and  £2.  For  example,  when  Buffer  1  operates 
under  a  constant-impedance  mode,  BESSi  supplies  any  shortfall 
in  power  for  £1  when  the  network  voltage  is  depressed.  This 
will  decrease  the  buffer  DC-link  voltage  udCi-  To  ensure  £1 
voltage  quality  is  not  affected,  however,  udCi  has  to  be  kept 
within  certain  limits,  denoted  by  the  range  [r>DC, lower,  ^dc, upper], 
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or  the  inverter  associated  with  Si  will  operate  outside  its  control 
range.  If  the  latter  case  occurs,  the  inverter  protection  system 
will  be  triggered  and  the  supply  to  Si  will  be  disrupted.  This 
situation  must  be  avoided.  The  extent  fdCi  will  deviate  from  its 
nominal  value  depends  very  much  on  the  capacity  of  BESSi. 
Similar  statements  can  be  made  concerning  BESS2.  Hence,  it 
is  necessary  to  determine  the  respective  BESS  capacity  in  order 
that  the  buffers  help  Si  and  S2  loads  ride  through  the  worst  fault 
conditions  of  the  scheme.  In  this  way,  the  power  buffers  will 
ensure  both  stability  enhancements  while  the  supply  quality  to 
the  loads  will  be  guaranteed.  Next,  a  method  to  determine  the 
BESS  capacity  will  be  described. 

2.4.  BESS  model  description 

Several  models  that  can  be  used  to  describe  the  behaviour 
of  a  battery  cell  have  been  proposed  by  researchers  [10-12]. 
The  model  proposed  in  reference  [10]  pertaining  to  a  lead-acid 
battery  cell  has  been  selected  for  incorporation  in  the  present 
buffer  scheme.  The  reason  for  the  selection  is  because  the  focus 
is  on  transient  stability  enhancement;  the  timescale  of  interest 
is  well  within  the  so-called  short-term  response  characteristics 
of  the  battery  cell  model  given  in  reference  [10].  The  adopted 
cell  circuit  model  is  reproduced  in  Fig.  2(a).  As  shown,  Ebxj  is 
the  cell  open-circuit  voltage,  internal  resistance  Rsjj  represents 
the  ohmic  resistance  due  to  polarization  in  the  cell,  and  Rpjj 
and  Cpxj  are  used  to  describe  the  discharge  transfer  loss  and 
dynamics  between  parallel  plates  of  electrodes  and  electrolyte, 
respectively. 

In  Fig.  2(b),  the  typical  response  of  the  cell  terminal  voltage 
Vb,i,j  to  a  step  cell  current  change  ibxj  is  shown.  It  is  noted  that 
a  step  increase  in  i^ij  would  result  in  the  cell  terminal  voltage 
Vbxj  decreasing  almost  instantaneously  before  it  approaches  its 
steady- state  value  in  an  exponential  manner.  In  describing  this 
behaviour,  the  change  in  battery  cell  terminal  voltage  due  to  the 
step  current  has  been  described  in  reference  [10]  in  the  following 
way: 

vb,i,j  =  Eb,i,j  —  Rs,i,jlb,i,j  ~  Rp,i,jib,i,j(  1  —  e  ^  ?’7)  (1) 

where  r/j  is  the  time  constant  of  battery  cell  and 

%j  =  RpfjCpfJ  (2) 

Typically,  r/j  is  of  the  order  of  0.2  s  to  a  few  seconds,  depend¬ 
ing  on  the  magnitude  of  the  step  current  change.  Table  1  shows 


Table  1 

Battery  parameters  used  in  simulations 


A/ (A) 

RS  (fl) 

Rp(Q) 

Cp(  F) 

t(s) 

25 

0.005 

0.015 

20 

0.3 

50 

0.008 

0.003 

900 

2.7 

100 

0.0052 

0.003 

1000 

3 

200 

0.0032 

0.002 

1000 

2 

400 

0.0022 

0.001 

750 

0.75 

600 

0.002 

0.0008 

600 

0.48 

800 

0.0018 

0.00075 

400 

0.3 

1000 

0.0015 

0.0007 

200 

0.14 

Fig.  2.  (a)  Schematic  showing  interconnection  of  lead-acid  battery  cells  in  BESS 
and  (b)  typical  response  of  battery  cell  terminal  voltage  Vb,ij  to  step  change  in 
cell  current  ib,ij. 


the  typical  parametric  values  of  the  cell  resistances  and  capaci¬ 
tance  [10]  and  the  corresponding  tq. 

Assume  the  BESS  has  m  cells  connected  in  series  in  each 
string  and  a  total  of  n  parallel  strings,  as  shown  in  Fig.  2(a). 
Assume  the  mn  cells  are  identical;  the  characteristic  of  the  ith 
series  cell  on  the y'th  parallel  string  of  the  BESS  is  described  by 
Eqs.  (1)  and  (2).  The  BESS  terminal  voltage  is  denoted  by  fdc 
since  the  BESS  is  connected  directly  across  the  DC-link  of  the 
buffer  system.  The  BESS  equivalent  circuit,  based  on  Fig.  2(a), 
is  given  in  Fig.  3,  where  the  BESS  internal  parameters  Eb,  Rs, 
Rp  and  Cp  can  be  expressed  in  terms  of  the  cell  parameters  as 
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1  • 

Fig.  3.  Equivalent  circuit  of  BESS. 


follows: 

Eb  =  mEb,ij/n  ' 

Rs  =  in  Rs  j  j/n 

”J/  ^  (3) 

Rp  —  m  Rpij  /  n 

Cp  =  nCpXj/m  , 

In  the  manner  similar  to  Eq.  (1),  the  BESS  terminal  voltage 
fdc(0  can  be  described  in  terms  of  the  BESS  parameters  and 
following  a  step  change  of  the  BESS  current  ib  as 

fdc  =  Eb  -  Rsib  Rpib(  1  e_^T)  (4) 


2.5.  DC -link  voltage  determination 


As  explained  in  Section  2.3,  it  is  important  that  the  DC-link 
voltages  must  be  maintained  within  limits.  Using  the  equiva¬ 
lent  circuit  model  developed  earlier,  the  BESS  terminal  voltage 
can  be  determined  as  follows.  In  the  study  of  transient  stability 
and  in  a  step-by-step  manner,  the  quantities  P\i{t),  Pgi(0  and 
Pg2(0  are  readily  calculated.  At  Buses  1  and  2,  Psl  and  Ps2  are 
known  constants,  whereas  Pdx  and  Po2  vary  with  the  square  of 
the  magnitude  of  the  respective  bus  voltage.  The  BESS  output 
powers  Pbx  ( t )  and  Pb2(t)  can  therefore  be  calculated  as  follows: 


Pbi(t)  =  PSl  ~  Pgi  +  P\2  +  Pdi  1 
Pb2(t)  =  PS2  -  Pg2  ~  +  Pd2  J 


Consider  one  of  the  buffers  and  omit  the  subscript  1  or  2  to 
simplify  the  expressions  that  follow.  Approximate  the  contin¬ 
uous  function  Pb(t)  over  the  interval  of  interest  T  by  a  series 
of  N+  1  values:  Pb( 0),  Pb( l),  •  •  •,  Pb(k),  •  •  •,  Pb(N)  at  times  f(0), 
t( i), . . .  t(k),  . . t(N).  Fig.  4  shows  the  ( k  —  l)th,  ^th  and  (k+  l)th 
intervals.  Assume  a  uniform  time  interval  so  that  each  segment 
is  of  duration: 


In  this  way,  Pbit)  may  then  be  approximated  by  N  pulses, 
each  of  duration  At,  and  the  magnitude  of  each  pulse  is  given 
by: 


Pb(k)  = 


(Pbjk- 1)  +  Pbjlc)) 
2 


k  =  1,2, N 


(7) 


Hence,  Pb(k)  will  be  the  power  pulse  applied  over  the  Mi 
time  interval.  Accordingly,  the  battery  current  step  change  due 
to  Pm  is: 


lb(k)  =  —  k  =  1,2, N 

VDC(k-l ) 


(8) 


where  r  is  the  time  constant  of  the  BESS  and 
r  =  RpCp  (5) 


where  fdc(&-1)  is  the  BESS  terminal  voltage  at  the  end  of  the 
( k  —  l)th  time  interval. 

Next,  r>DC (k)  can  be  generated  through  the  following  step-by- 
step  and  iterative  method.  Consider  udc  at  the  end  of  the  first 


The  circuit  model  of  the  BESS  is  amendable  for  analysis, 
except  for  one  complication.  In  adopting  the  equivalent  circuit 
model,  it  has  been  stated  earlier  that  the  cell  resistance  and  capac¬ 
itance  values  are  not  constant  but  are  functions  of  the  magnitude 
of  the  current  step.  In  fact,  the  state-of-discharge  of  the  cell  and 
temperature  at  which  the  cells  operate  at  also  affect  the  paramet¬ 
ric  values.  Fortunately,  the  latter  two  factors  could  be  considered 
constant  over  a  fault  disturbance  event  because  the  duration  of 
the  fault  and  its  clearance  is  usually  too  short  for  the  battery 
temperature  and  discharge  state  to  vary  appreciably.  In  this  way, 
the  values  of  the  parameters  would  only  depend  on  the  magni¬ 
tude  of  the  discharge  current  step.  In  what  follows,  it  will  be 
assumed  that  the  parametric  values  of  Rs ,  Rp  and  Cp  would  be 
known  for  given  magnitude  of  the  current  step.  The  information 
can  be  obtained  either  from  the  cell  manufacturer  or  through 
experimental  measurements. 
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time  step,  i.e.,  k  =  1 .  The  initial  BESS  terminal  voltage  is  given  by 
^DC(O)  =  Eb,  assuming  the  battery  current  is  zero  initially  and 
the  BESS  open-circuit  voltage  is  known.  As  shown  in  Fig.  5, 
assume  the  first  power  pulse,  Pb(i),  is  positive.  From  Eq.  (8),  the 
corresponding  battery  current  step  change  is  given  by: 

?«,,  =  ^  (9) 

H)C(O) 

With  the  known  step  current,  one  can  utilize  the  correspond¬ 
ing  BESS  parametric  values,  denoted  as  Rs(i),  Rp( l)  and  Cp( i),  to 
determine  the  DC -link  voltage  at  the  end  of  the  first  time  interval 
using  Eq.  (4),  i.e., 

^DC(i)  =  *T>C(0)  -  ^(1)4(1)  ~  ^(1)4(1)(1  -  e-(A*/T(1)))  (10) 


where  r(i)  =  Rp(\)Cp(iy 

The  net  fdc  variation  Ai;dc(1)  at  the  end  of  the  first  time  step 
is  derived  from  Eq.  (10),  i.e., 

Al^DC(l)  =  ^DC(l)  -  ^DC(O)  =  ^VDC(1,1) 

=  -(^(1)4(1)  +  rp( i)4(i)(l  -  e_(A^T(1))))  (1 1) 

where  Av^c^  ^  denotes  the  voltage  variation  caused  by  the  first 
positive  current  pulse  at  the  end  of  the  first  time  interval  f  i). 

Fig.  5  shows  that  r>DC  varies  over  the  first  time  interval.  Since 
Pb( l)  is  assumed  constant  over  the  interval  At,  the  battery  current 
will  in  turn  vary  with  udc-  Considering  this  factor,  an  inner-loop 
calculation  step  is  suggested  in  the  numerical  solution  method 
step  to  improve  the  accuracy  of  the  DC-link  voltage  calcula¬ 
tion.  The  inner-loop  calculation  is  described  as  follows.  Using 
AfDC(i)  obtained  from  Eq.  (10)  as  an  intermediate  variable,  one 
can  derive  the  corresponding  battery  current  as: 


^(l) 


vdco) 


(12) 


Then  update  lb  by  taking  the  average  value  of  4  from  (9)  and 
Th{V)  fromEq.  (12).  One  then  recalculates  Ar>DC(l)  using  Eq.  (11) 
and  fDC(l)  through  Eq.  (10)  with  the  new  lb .  The  above  improved 
battery  current  and  terminal  voltage  calculation  procedure  is 
repeated  until  udc(1)  convergences.  Since  it  has  been  known 


that  the  battery  time-constant  r  is  much  larger  than  the  step- size 
At,  it  is  envisaged  that  this  inner-loop  calculation  will  only  result 
in  relatively  minor  corrections  in  lb  during  the  iterative  process. 
Hence,  convergence  can  be  attained  rather  quickly. 

From  the  start  of  the  second  time  interval,  i.e.,  t>t( i),  Pb(2) 
is  introduced  as  shown  in  Fig.  6(b).  A  new  battery  current  step 
change  lb(2)  is  introduced  where 


4(2)  = 


VDC(l) 


H)C(l)  is  obtained  from  Eq.  (10).  To  simulate  completely  the 
net  effect  of  the  first  pulse,  however,  it  is  necessary  to  combine 
the  effects  shown  in  Fig.  6(a)  with  those  of  Fig.  6(c),  i.e.,  to 
introduce  a  negative  current  step  of  magnitude  lb  at  the  start 
of  the  second  time  interval,  as  shown  in  Fig.  6(c).  This  negative 
step  current  produces  a  voltage  response  at  the  end  of  the  second 
time  interval  thus 


aUxxi,2)  —  Rs(  1)4(1)  +  ^/?(i)4(i)(l  -  e  (Ar/Td)))  (i3) 


where  Av^c^  2)  denotes  the  voltage  variation  caused  by  the 
negative  step  current  at  the  end  of  the  second  time  interval  t( 2). 
In  Eq.  (13),  it  is  assumed  the  negative  step  current  has  the  same 
BESS  parametric  values  as  that  of  a  positive  step  current  of  the 
same  magnitude. 

Hence,  at  the  end  of  the  second  time  step,  the  BESS  voltage 
variations  caused  by  the  positive  step  currents  7 b  and  lb  are: 

Avdc(i,2)  =  -(Rsmibm  +  flp(i)4(i)(i  -  e-(2Ar4a>>)) 


A"DC(2,2)  =  -(^(2)4(2)  +  tfp(2)4(2)(l  ~ 

where  Rs(2),  Rp( 2)  and  Cp(2)  are  the  BESS  parameters  cor¬ 
responding  to  a  current  step  change  of  magnitude  lb  and 
T(2)  =  Rp(2)Cp(2).  At  the  same  time,  the  negative  step  current  pro¬ 
duces  the  voltage  change  Av^c^  2\  as  shown  in  Eq.  (13).  Hence, 
the  battery  terminal  voltage  variation  at  the  end  of  the  second 
time  step  is: 

Ar>DC(2)  =  Ai;qC(;i?2)  +  Av^c(2,2)  +  ^UDC(1,2) 
and  r>DC  at  f  2)  is 


Pb(D  . . 

t(0)  fl)  t 

Vdc  ^ 

Vdc(0)  - 1 . J  -  , 

L,  f  A4/cm  =Av<fc(i>i) 

o' - * - i - ► 

t(o)  fi)  t 

Fig.  5.  Estimation  of  DC-link  voltage  at  end  of  first  time  interval. 


FDC(2)  =  UXXO)  +  AfDC(2)  (14) 

Having  obtained  fdc(2)  from  Eq.  (14),  the  same  inner-loop 
updating  procedure  described  earlier  for  7^  to  calculate  A^dc 
over  the  first  time  interval  can  also  be  applied  until  fdc(2)  con¬ 
vergences. 

The  above  fdc  calculation  procedure  is  repeated  over  the 
interval  T.  In  general,  at  the  kth  interval,  vr>C(k)  is  determined  by 
the  following  equations: 

k  k- 1 

VDC(k)  =  VDC(0)  +  AuDC (M,k)  (15) 

M=  1  M=  1 

where 

AvDC  (M,k)  =  ~iRS(M)h(M)  +  Rp(M)h(M) 
x  (t  _  Q-((k-M+l)At/r (M)))) 


(16) 
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Fig.  6.  Estimation  of  DC-link  voltage  at  end  of  second  time  interval. 


Ai>dc (M,k)  =  RS(M)h(M)  +  Rp(M)Ib(M)0  -  e-^-^'/W) 

(17) 

Rs(m ),  Rp(M)  and  CP(M)  are  the  BESS  parameters  correspond¬ 
ing  to  a  current  step  ~h(M)  and  T(m)  =  Rp(M)Cp(m >  k  =  2,  3,  . . 
N. 

Combining  Eqs.  (10)  and  (15),  one  can  obtain  general  expres¬ 
sions  for  determining  i>dc  at  the  end  of  the  &th  time  step  as: 

"  udc(0)  +  Addc(1)  k=  1 

k  k—  1 

^DC(&)  —  V — >  V — >  _  (18) 

«DC(0>  +  2_^AVDC(M,k)  +  2_^AVDC(M,k)  k  =  2,3...N 

,  M=  1  Af=  1 

where  Ai;dc(1)?  ^^dc(m  k)  anc*  ^^dccm^)  are  calculated  using 
Eqs.  (11),  (16)  and  (17),  respectively. 

2.6.  BESS  capacity  determination 

It  is  now  clear  that  the  BESS  not  only  compensates  for  the 
balance  of  the  power  between  the  constant-power  load  demand 
and  that  supplied  by  the  upstream  utility  system,  but  also  ensures 
that  supply  quality  to  the  load  does  not  deteriorate. 

Also,  as  explained  in  Section  2.2,  when  a  fault  occurs  and 
the  buffer  is  switched  to  the  constant-impedance  control  mode, 


power  has  to  be  extracted  from  the  buffer  BESS  to  meet  the  bal¬ 
ance  of  the  power  demanded  by  the  downstream  constant-power 
load.  Discharging  power  from  the  BESS  will  consequently 
depress  the  BESS  terminal  voltage  fdc-  If  the  BESS  capacity 
is  insufficient,  fdc  may  be  forced  to  drop  below  its  set  lower 
limit  fdc, lower-  As  described  in  Section  2.3,  this  is  not  accept¬ 
able.  One  could  increase  the  BESS  capacity  to  prevent  fdc 
from  falling  below  fdc, lower-  This  can  be  achieved  by  adding 
additional  strings  of  battery  cells.  From  Section  2.4,  increas¬ 
ing  the  number  of  strings  n  will  correspondingly  reduce  the 
internal  resistances  Rs  and  Rp ,  whereas  the  capacitance  Cp 
will  increase.  For  a  given  P &,  the  discharge  current  will  be 
reduced,  and  hence,  the  terminal  voltage  deviations  will  be 
less. 

Over  the  transient  interval  when  the  power  system  read¬ 
justs  the  power  flows  and  network  voltages,  it  is  also  likely 
that  the  amount  of  power  supplied  by  the  upstream  source  can 
be  greater  than  that  demanded  by  the  downstream  constant- 
power  load.  Under  such  circumstances,  the  surplus  power  will 
be  absorbed  by  the  BESS  and  fdc  increases.  Hence,  the  capacity 
of  BESS  should  also  be  sufficiently  large  to  guarantee  fdc  does 
not  exceed  the  pre-set  upper  limit  fdc, upper-  Again,  this  can  be 
achieved  by  increasing  the  number  of  parallel  branches  of  the 
cells.  The  capacity  of  the  BESS  may  therefore  be  determined 
based  on  the  worst  fault  disturbance  condition  for  which  the 
buffer  scheme  is  designed. 
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3.  Simulation  results 

To  illustrate  the  design  procedure  of  the  previous  section,  the 
simple  power  system  depicted  in  Fig.  1  was  simulated  using  the 
well-known  package  PSCAD.  In  this  example,  a  single-circuit 
radial  link  between  Buses  1  and  2  is  assumed.  The  line  and  gen¬ 
erator  parameters  used  for  the  simulations  are  listed  in  Appendix 
(A).  Data  is  taken  from  references  [13,14].  The  AC1A  type 
excitation  and  mechanical-hydraulic  controlled  (MHC)  gover¬ 
nor  systems  [14]  are  also  included  in  the  generator  systems. 
The  power  flow  pattern  at  a  nominal  condition  is  such  that  G\ 
has  a  generation  of  1575  MW  and  G2  is  of  75  MW.  The  load 
levels  at  Buses  1  and  2  are  1500  and  150  MW,  respectively. 
Hence,  75  MW  is  imported  from  Bus  1  into  Bus  2  through  the 
transmission  link. 

Assume  half  of  the  loads  at  Buses  1  and  2  are  essential  (super¬ 
premium)  and  are  to  be  protected  by  the  power  buffers.  By 
assuming  a  line  fault  condition,  the  response  of  the  power  sys¬ 
tem  was  studied  and  from  the  power  flows  and  Eq.  (6),  P &  on  the 
BESS  can  be  obtained  for  both  the  buffers.  fdc  of  the  respective 
buffer  is  calculated  using  the  computational  method  described 
in  Section  2. 

3.1.  A  temporary  mid-line  three-phase  fault 

A  temporary  10-cycle  three-phase-to-ground  fault  event  at 
the  midpoint  of  the  transmission  line  is  examined  first.  By 


varying  the  buffer  switching  pattern  to  realize  either  constant- 
impedance  or  constant-power  modes  on  the  basis  of  the 
variations  of  P\2  about  the  pre-fault  transfer  level  of  75  MW,  it 
is  intended  to  demonstrate  that  the  stability  of  the  power  system 
is  improved.  As  stated  earlier,  Power  Buffer  2  is  switched  into  a 
constant-impedance  mode  during  the  fault  occurrence.  BESS2 
discharges  to  compensate  for  the  mismatched  power  between  the 
constant  power  S2  demanded  and  that  supplied  by  the  upstream 
sources.  The  power  profile  of  BESS2  is  calculated  using  Eq.  (6) 
and  is  given  in  Fig.  7(a),  in  which  positive  values  denote  power 
discharging  from  BESS2  while  negative  ones  indicate  BESS2 
is  being  charged.  At  Buffer  1,  there  is  a  negligible  amount  of 
power  exchange  between  BESSi  and  the  external  system.  This 
is  because  the  constant  mechanical  power  generated  by  G\  is 
greater  than  that  demanded  by  the  load  connected  to  Bus  1. 
Hence,  Buffer  1  operates  under  a  constant-power  mode  during 
the  fault  to  reduce  acceleration  of  the  G\  rotor.  Therefore,  only 
the  performance  and  simulation  results  of  Buffer  2  are  described 
as  follows. 

The  BESS  cells  have  the  parametric  values  given  in  Table  1. 
Assume  the  modulation  indices  of  the  power  buffers  are  at  the 
typical  value  of  0.8  and  the  initial  DC-link  voltage  of  Power 
Buffer  2  is  269.44  kV.  Therefore,  each  BESS  branch  will  be 
composed  of  40,825  cells  since  each  cell  is  rated  at  6.6  V.  Using 
the  computational  procedure  described  in  Section  2.5  and  based 
on  Pb2  shown  in  Fig.  7(a),  fdc  is  calculated.  The  time  step  used 
is  the  half-cycle  of  the  power  frequency.  The  r>DC  so  obtained  is 


Fig.  7.  BESS2  response  following  temporary  10-cycle  three-phase  mid-line  fault:  (a)  P&2,  (b)  i>dc2  and  (c)  rotor-angle  difference  (6\  —  Of)  between  G\  and  G2:  with 
and  without  power  buffers. 


X.Y.  Wang  et  al.  /  Journal  of  Power  Sources  179  (2008)  819-829 


827 


shown  in  Fig.  7(b).  In  these  plots,  the  base  power  and  DC-link 
voltage  values  used  are  100  MW  and  269.44  kV,  respectively. 

For  satisfactory  operations  of  the  downstream  protected 
loads,  assume  the  DC-link  voltage  is  constrained  to  operate 
within  the  range  [0.9,  1.05]  p.u.  This  is  a  reasonable  assump¬ 
tion;  for  example,  adjustable- speed  drive  systems  would  only 
operate  satisfactorily  within  this  range  [15].  Assume  BESS2  has 
only  one  string  of  the  cells.  From  Fig.  7(b),  it  is  seen  that  in 
supporting  the  load  demand  S2,  ^dc2  is  depressed  below  the  set 
lower  limit  over  the  initial  stage  of  the  fault  disturbance.  In  this 
case,  S2  will  not  be  supplied  at  a  voltage  level  that  is  acceptable 
for  satisfactory  continuous  operation. 

One  can  attempt  to  increase  BESS2  capacity  by  adding  one 
more  string  of  the  cells.  The  corresponding  t>DC2  Pr°file  f°r  the 
two-string  BESS2  is  also  shown  in  Fig.  7(b).  The  two-branch 
BESS2  is  observed  to  be  capable  of  maintaining  t>DC2  well  within 
the  limits.  Therefore,  the  capacity  of  BESS2  is  determined  to 
consist  of  two  strings  because  in  this  instance  the  DC-link  volt¬ 
age  is  improved  and  the  voltage  quality  of  S2  is  of  an  acceptable 
level. 

The  effect  of  the  power  buffers  in  enhancement  of  the  tran¬ 
sient  stability  can  be  seen  from  Fig.  7(c).  As  shown  by  the  dotted 
line  that  corresponds  to  the  case  without  buffers  incorporated, 
the  rotor-angle  difference  continues  to  increase  monotonically 
and  the  system  loses  synchronism  consequently.  When  the  pro¬ 
posed  buffer  operation  is  utilized,  however,  instability  is  avoided 
as  can  be  seen  from  the  solid  line  in  Fig.  7(c).  It  clearly  demon¬ 


strates  that  the  buffer  systems  do  function  as  effective  network 
system  dampers  or  stabilizers. 

An  interesting  observation  can  be  made  from  Fig.  7(b).  It  is 
seen  that  the  udc  post-fault  steady-state  level  remains  the  same 
as  that  of  the  pre-fault  level,  although  it  is  clear  from  Fig.  7(a) 
that  there  is  a  net  export  of  energy  from  BESS 2  to  the  load  over 
the  period  of  study.  With  the  net  discharge  of  energy  from  the 
BESS,  the  battery  terminal  voltage  is  expected  to  decrease  in 
the  longer  term  because  its  internal  emf  (. E &)  decreases  as  the 
BESS  state-of-discharge  (SoD)  increases  [10].  In  the  present 
work,  it  has  been  assumed  that  E b  remains  constant  over  the  fault 
incident.  For  longer  term  study  of  (say)  several  tens  of  seconds, 
it  is  necessary  to  include  the  SoD  effect  into  the  analysis. 

3.2.  A  temporary  three-phase  fault  close  to  Bus  2 

The  capacity  of  the  BESS2  determined  for  the  midpoint  fault 
condition  is  evaluated  for  its  suitability  for  the  same  temporary 
three-phase  10-cycle  fault  event,  but  is  now  located  at  a  distance 
corresponding  to  1  km  from  Bus  2.  Pb2  for  this  case  is  given  in 
Fig.  8(a),  while  the  noc  computed  is  shown  in  Fig.  8(b). 

Comparing  P^  shown  in  Fig.  8(a)  with  that  in  Fig.  7(a),  it  is 
observed  that  with  the  fault  occurring  closer  to  Bus  2,  BESS2  has 
to  increase  its  discharged  power.  Indeed,  as  shown  in  Fig.  8(b) 
and  under  such  a  fault  condition,  BESS2  with  one  cell  string 
would  result  in  i>dc2  not  only  falling  below  the  lower  limit  over 
the  fault  interval  but  also  subsequently  rising  momentarily  above 


time(s) 


Fig.  8.  BESS2  response  following  temporary  10-cycle  fault  1  km  from  Bus  2:  (a)  Pb2,  (b)  udc2  and  (c)  rotor-angle  difference. 
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Fig.  9.  BESS2  response  of  100%  of  constant  power  loads  connected  at  both  ends  under  temporary  three-phase  fault  1  km  from  Bus  2:  (a)  Pt2,  (b)  l>dc2  and  (c) 
rotor- angle  difference. 


its  upper  limit  during  the  system  recovery  period.  The  estimated 
fdc2  when  the  BESS2  consists  of  two  cell  strings  is  also  given 
in  Fig.  8(b)  and  is  seen  to  produce  satisfactory  performance. 
Again,  as  shown  in  Fig.  8(c),  the  power  system  is  transiently 
stable  with  such  a  load  buffer  scheme. 

3.3.  Further  discussion  of  results 

One  more  interesting  issue  pertains  to  the  proportion  of  the 
super-premium  loads.  Instead  of  assuming  the  load  connected 
to  each  bus  should  consist  of  equal  amounts  of  constant-power 
super-premium  load  and  constant-impedance  load,  additional 
studies  have  been  carried  out  when  100%  of  the  loads  are  taken 
to  be  super-premium  and  are  to  be  protected  by  the  buffers. 
Hence,  the  pre-fault  load  capacity  of  S\  is  1500  MW  and  that  of 
S2  is  150  MW. 

For  the  purpose  of  comparison,  the  simulation  studies  corre¬ 
sponding  to  the  above-described  load  pattern  under  temporary 
10-cycle  three-phase  fault  cases  at  mid-line  and  1  km  from  Bus  2, 
respectively,  are  repeated.  Due  to  the  limitation  of  space,  only  the 
results  for  a  fault  that  is  1  km  from  Bus  2  are  given  in  Fig.  9.  For 
this  fault  case,  the  power  profile  of  BESS2  is  shown  in  Fig.  9(a), 
from  which  it  is  found  that  the  output  power  of  BESS2  has  to  be 
increased  due  to  the  increase  in  the  proportion  of  the  constant- 
power  load  at  Bus  2.  Therefore,  the  capacity  of  BESS2  has  to  be 
enlarged.  This  point  can  be  seen  from  the  DC-link  profile  given 


in  Fig.  9(b),  where  even  for  the  two-string  BESS,  fdc  still  vio¬ 
lates  the  lower  and  upper  udc  limits  during  both  fault  and  post 
fault.  One  additional  string  is  added  to  BESS2  and  the  resulting 
r>DC  profile  is  shown  as  the  solid  line  in  Fig.  9(b).  The  improved 
transient  stability  with  the  help  of  the  power  buffer  can  be  seen 
from  the  rotor- angle  difference  profile  in  Fig.  9(c). 

From  the  results  of  the  numerical  examples,  it  is  not  surpris¬ 
ing  to  note  that  as  the  fault  event  is  moved  electrically  closer 
to  the  load-buffer  or  when  the  proportion  of  the  super-premium 
loads  is  increased,  it  is  necessary  to  increase  the  BESS  capacity 
to  cater  for  the  constant-power  demand  of  the  loads.  There¬ 
fore,  the  capacities  of  the  BESSs  within  the  power  buffers  have 
to  increase  with  increasing  penetration  of  the  super-premium 
constant-power  loads  in  the  network. 

4.  Conclusions 

From  the  basic  function  of  power  buffers  for  voltage  quality 
improvement,  an  extension  of  the  buffer  application  to  enhance 
power  system  stability  has  been  explored.  The  focus  of  the  inves¬ 
tigation  is  on  the  design  of  the  BESS  of  the  buffers.  When  a 
buffer  operates  under  a  constant-impedance  mode,  the  BESS  dis¬ 
charges  or  absorbs  power.  This  will  cause  its  terminal  voltage  to 
deviate.  A  computational  technique  to  determine  the  terminal 
voltage  has  been  developed  based  on  the  battery  power  dis¬ 
charged  or  absorbed.  As  the  BESS  is  connected  across  the  buffer 
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DC-link,  the  voltage  deviations  will  have  an  impact  on  the  volt¬ 
age  quality  of  the  downstream  protected  load.  The  proposed 
method  allows  a  suitable  capacity  of  the  BESS  to  be  determined 
so  that  acceptable  voltage  quality  for  the  protected  load  will  be 
obtained.  Numerical  results  show  that  the  capacity  of  the  BESS 
increases  with  the  amount  of  the  protected  load  and  the  severity 
of  the  fault.  Through  the  implementation  of  the  proposed  power 
buffer  scheme,  the  stability  of  the  network  will  also  be  enhanced 
under  system  disturbance  conditions. 

Appendix  A 

A.l.  (A)  Line  parameters  used  in  simulations  [13] 

System  voltage  =  132  kV,  nominal  frequency  =  50  Hz, 
line  length  =  200  km,  r  =  0.068  £2  km-1,  v  =  0.404  Q  km-1, 
Z?  =  7.59E— 6Skm-1. 

A.  2.  (B)  Equivalent  machine  parameters  used  in 
simulations  [14] 

Gi:  SN  =  2000MVA,  //  =  4.0s,  l.OOp.u.,  X'd  = 

0.30p.u.,  X"d  =  0.20p.u.,  2Q  =  0.60p.u.,  X"q  =  0.20p.u., 
Xa  =  0.12p.u.,  r'0  =  5.00s,  r"0  =  0.08s,  T'0  =  -, 

**,=  0.12  s. 

G2:  Sn  =  100  MVA,H= 4.0  s,Xd  =  2.00 p.u.,  X'd  =  0.25  p.u., 
X"d  =  0.20  p.u.,  X?  =  0.50  p.u.,  X"q  =  0.20  p.u.,  Xa  =  0.15p.u„ 
Td0  =  6.00  s,  Td0  =  0.03  s,  7^0  =  0.60  s,  r"0  =  0.06  s. 

A. 3.  ( C)  AC1A  type  exciter  and  regulator  parameters  used 
in  simulations  (p.u.)  [14] 

Tr  =  0.0,  Ka  =  400.0,  Ta  =  0.02,  TB  =  0.0,  Tc  =  0.0,  KF  =  0.03, 
Tf=  1.0,  Ke=  1.0,  7£  =  0.8,  ATo  =  0.38,  Kc  =  0.2,  V«MAX  = 
-6.03,  V£min  =  -5.43,  VAmax  =  14.5,  VAmin  =  -14.5. 


A.  4.  (D)  MHC  governor  parameters  used  in  simulations 
[14] 

Gi:  TP  =  0.05  s,  ^  =  5.0,  TG  =  0.2  s,  7?P  =  0.04, 

Rt  =  0A ,  Tr  =  5.0,  Gmax  =  1.0  p.u.,  Gmin  =  0.0  p.u., 

^max, open  —  0.16  p.U.  S  7?max  ciose  =  0.16  p.U.  S 

G2:  Kg  =  20,  Tsr  =  0As,  Tsm  =  0.2s,  TSi  =  0.2s,  Lc i=0.2, 
Lc2  =  -0.5,  L/i  =  0.2,  Ln  =  -0.5. 
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